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OPTIMAL REUSABLE-TUG AND EXPENDABLE-KICKSTAGE TRAJECTORIES FOR HIGH- 
ENERGY PLANETARY MISSIONS INCLUDING CORRECTION FOR NODAL PRECESSION 

by Janos Borsody 
Lewis Research Center 

SUMMARY 

Equations are derived by using branched trajectx)ry optimization techniques and the Maximum 
Principle to design trajectories that maximize the payload capability of a reusable, Space Shuttle 
upper stage (Tug) /expendable kickstage configuration for planetary missions. The analysis takes 
into account the Space Shuttle (SS) orbit nodal precession. The Tug/kickstage trajectories are de- 
signed to correct for this nodal precession. Since the kickstage is expended after payload injection, 
ill nodal correction maneuvers are done by the Tug. Planetary injection conditions that must be 
satisfied at kickstage burnout are defined in terms of the mission energy and the declination and 
eight ascension of the outgoing asymptote. 

The multipoint boundary value problem resulting from the mathematical analysis has been 
solved by using a Newton-Raphson iteration technique. Partial derivatives of the boundary condi- 
.ions required by this technique are approximated by perturbing the initial conditions one at a time, 
ntegrating the equations of motion and the adjoint variational equations, and observing the changes 
n boundary conditions. The amount by which each initial condition is perturbed is arbitrary; how- 
3ver, in obtaining the sample case results, it was found that a 0. 1 percent perturbation gave suffi- 
ciently accurate partial derivatives. Using these partial derivatives and the Newton-Raphson iter- 
itor led to a stable and uniform convergence. 

Numerical results are presented for a reusable cryogenic Tug design based on current state-of- 
.he-art technology and an expendable, space- storable, pressure- fed kickstage. Detailed analyses 
)f these trajectories are presented for two mission energies. In addition, the effects of several 
nission and stage parameters on payload capability are also investigated. 

Nominal mission characteristics used for the sample cases analyzed include a 185-kilometer 
circular departure and rendezvous orbit with an orbital inclination of 28. 5®. The Tug is constrained 
to return to the Space Shuttle (SS) within 1 day. It is further assumed that the Tug corrects for all 
nodal precession of the SS orbit during the outbound burn and the retrobum. 

A 10-minute coast time after the first Tug burn is included to allow for safe separation of the 
expendable kickstage and payload from the Tug. At the end of this coast phase, the kickstage ignites 
and accelerates the payload to the desired injection conditions. 

Detailed performance analysis was done for two energies, one typical of a comet mission 
(50 km^/sec^) and the other representing an outer planet mission (100 km^/sec^, for declinations of 
the outgoing asymptote between -34® and +34®. Ideal impulsive performance calculations are shown 
to be highly optimistic for this mission because large gravitational and thrust-to-weight losses are 
encountered during the flight. Payload loss to correct for nodal precession is snfiall for both missions. 

A parametric analysis of several stage and mission parameters was also performed for the outer 
planet mission. The major mission and Tug parameters that were perturbed include initial SS orbital 
inclination, total trip time, Tug/kickstage separation coast time, and Tug and kickstage engine thrust 
levels. These results show that payload is very sensitive to thrust level and kickstage separation 
coast time. Optimum total Tug trip time for this mission is approximately 26 hours; however, the 
payload penalty for trip times between 20 and 36 hours is quite small 

For a declination range between -34® and +34®, maximum performance is obtained for Tug 
departure-orbit inclinations of 28. 5°. Higher declinations can be reached by increasing the SS orbital 
inclination with small payload penalty. 



INTRODUCTION 


At the present time the National Aeronautics and Space Administration (NASA) is 
developing a reusable Space Shuttle (SS). The SS will deliver many different payloads 
to low Earth orbit. For missions beyond SS capability, a propulsive stage (or stages) 
along with the payload will be delivered and deployed by the SS into a low Earth orbit, 
and the propulsive stage (or stages) will inject the payload to its final orbit. This pro- 
pulsive stage configuration can be expendable, reusable, or a combination thereof. 

Previous NASA studies indicate that a substantial cost savings can be realized by 
using a reusable orbital propulsive stage. For the numerical example used in this re- 
port, a conceptual orbital propulsion stage has been selected, which will be referred to 
throughout this report as the ”Tv^. " The projected cost savings for the Tug are partly 
due to retrieval, refurbishment, and reuse of spacecraft and partly achieved throi^h the 
reuse of expensive propulsive-stage hardware. For mission energies beyond reusable 
Ti^ capability, either the expensive Tug hardware must be expended or multistage orbi- 
tal operations have to be considered. 

Performance of a single-stage reusable Tug for low-energy, inner planet missions 
has been analyzed in reference 1. In the present report a method is presented to deter- 
mine the maximum performance capability of a reusable Tug and an expendable kick- 
stage for medium-energy (comet) and high-energy (outer planet) missions. This con- 
figuration could be attractive since, instead of expending an expensive piece of hardware 
(the reusable Tug), only a less-expensive kickstage is expended to do the mission. It 
will be assumed that the Tug, expendable kickstage, and payload are launched with a 
single SS flight. This assumes that the Tug propellant load is varied with payload and 
kickstage weight to remain within the SS lift capability of 29 484 kilograms. 

Previous performance analyses of a reusable/expendable stage configuration have 
been largely limited to ideal impulsive computations, which assume infinite Tug and 
kickstage thrust. It has been further assumed in these studies that each burn occurs in 
low Earth orbit at an optimum point along the orbit. In reality, all burns have finite 
length (in some cases, several minutes), and some of the burns are performed far from 
Earth (depending on the trajectory constraints), which results in substantial gravita- 
tional losses. Some planar integrated performance data have been presented for multi- 
stage vehicle configurations in reference 2 without corrections for nodal precession. In 
the present report the effects of finite thrust level as well as nodal correction for the 
SS orbital precession are included in the analysis. 

Nodal precession is caused by the Earth’s oblateness and to a first-order approxi- 
mation it is a function of orbital inclination, eccentricity, semilatus rectum, and the 
time spent in orbit. Since the duration of this mission (the time from Tug departure to 
return to the SS orbit) is of the order of 1 day, the selected SS orbit processes approx- 
imately 8°. The Tug will have to correct for this nodal precession in order to return 
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to rendezvous with the SS. Nodal precessions of the Tug orbits are negligible because 
of the large semilatus rectum and eccentricity of these orbits. If nodal precession of 
the SS orbit is ignored as in reference 2, it can be shown that the Tug trajectories be- 
come planar and, furthermore, that the payload becomes independent of declination. 

This greatly simplifies the analysis and the solution of the two-point boundary value 
problem. 

The purpose of the present report is twofold. One is to derive the equations that 
govern this complex mathematical optimization problem; the other is to apply the re- 
sulting equations to a sample case and to evaluate the effects of nodal change on payload 
and flight profile. 

The mathematical equations are derived by using branched trajectory optimization 
techniques and the Maximum Principle. The multipoint boundary value problem result- 
ing from this analysis is solved by using a Newton- Raphson iterator. Partial deriva- 
tives of the boundary conditions with respect to the free initial conditions required by 
this iterator are obtained from finite difference perturbations. 

The method derived in this report was used to calculate the maximum payload 

capability for a fixed vehicle configuration and two mission energies of interest. The 
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first energy (50 km /sec ) is typical of a high-energy inner planet, comet, or asteroid 
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mission; the second (100 km /sec ) is a typical outer planet mission energy require- 
ment. 

The nominal vehicle selected for the sample case analysis consists of a cryogenic 
Tug and a space- storable kickstage. The Tug performance is based on a design that 
uses hydrogen and oxygen propellants; the pressure-fed, kickstage performance is 
based on a design that uses fluorinated oxygen and monomethyl hydrazine (MMH) pro- 
pellants. 

The effects of some trajectory parameters such as departure- orbit inclination, 
total trip time, and Tug/kickstage separation coast time on payload are investigated and 
discussed. The effect of Tug and kickstage thrust on payload is also investigated. 

These parametric studies were done only for the outer planet mission. 


ANALYSIS 

In the following sections the trajectory profiles are discussed in detail. Assump- 
tions and constraints associated with the problem are presented, and the mathematical 
optimization problem is formulated in the framework of the Maximum Principle. From 
this formulation and the transversality conditions, a set of auxiliary variational condi- 
tions are derived. These variational conditions are further reduced to give a set of 
boundary conditions that must be satisfied in order to maximize payload. A short dis- 
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cussion of the procedure used to solve the two-point boundary value problem resulting 
from the analysis is also included. 


Flight Profile and Simplifying Assumptions 
A planar trajectory profile is presented in sketch (a) for the mission to be analyzed. 


Perigee correction burn (AVg) 


Space shuttle V 
circular orbit— ' 


Circularization 
burn (AVp) 



/^Outgoing 
/' ' asymptote 


Expendable kickstage burn 


■/'-Outbound Tug burn 


(a) 


The flight profile consists of four burns by the Tug (outbound burn, retroburn, and 
apogee and perigee burns) and a single kickstage burn. The outbound Tug burn takes 
the kickstage and payload past Earth escape. This burn is followed by a fixed coast 
phase, which allows the Tug and kickstage to separate and the Tug to turn around to 
orient it for a retroburn. Retroburn is necessary because at separation the Tug is on an 
escape hyperbola and in order for it to return to the SS an elliptic return orbit must be 
established. Since during retroburn the Tug engine is pointed toward the kickstage and 
payload, plume impingement (causing possible payload contamination or damage to sen- 
sitive instrumentation) could be a problem. Also the highly corrosive exhaust products 
of the space -storable kickstage could damage the Tug, making refurbishment of this 
stage more difficult and costly. To alleviate this problem, a separation velocity must 
be imparted to the kickstage and payload, and the coast time must be kept long enough 
to allow for a safe separation distance to be achieved before either the Tug or the kick- 
stage are ignited. The required separation distance for this configuration has not been 
established and a fixed coast time is assumed for this analysis. At the end of this coast 
phase the Tug performs the retroburn, while the kickstage accelerates the payload to its 
given planetary injection conditions. Retroburn is terminated primarily to satisfy a 
given trip time. The Tug is also constrained to be in the same orbital plane at retro- 
burnout as the SS will be in at rendezvous. This constraint was introduced to correct 
for nodal precession of the SS orbit. 
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Nodal precession is caused by the Earth's oblateness, which induces a retrograde 
rotation of the SS orbit about the polar axis. This rotation is about 8 degrees per day 
for the selected SS orbit. Since the total Tug trip time (time from Tug departure from, 
to return to, the SS orbit) is about 1 day, the Tug must make an 8° plane change during 
the outbound burn and the retroburn. The flight following the retroburn will be planar. 

At the apogee of the intermediate orbit, established at retroburnout, a small perigee 
correction burn is performed followed by a perigee burn to accomplish the rendezvous 
with the SS. These burns are assumed to be impulsive and to occur exactly at apogee 
and perigee in the analysis. Gravitational and thrust- to- weight losses should be very 
small for these two burns since the burns are short and the thrust-to-weight ratio is 
quite high. 

The outbound burn, retroburn, and kickstage burn are numerically integrated, but 
the apogee and perigee burns are done impulsively. The assumption of impulsive apogee 
and perigee burns is fnade in order to simplify the mathematical optimization problem. 

It affects the problem in two ways: First, it allows these burns to be included analyti- 
cally, which reduces by four the number of variational boundary conditions that must be 
satisfied. Second, the equations of motion and the adjoint variational equations must be 
integrated only to the end of the retroburn (a small fraction of the total Tug trip time). 
This greatly reduces the effects of roundoff errors and error propagation and conse- 
quently improves the accuracy of the finite-difference partial derivatives required by 
the Newton-Raphson iteration. The flight during the integrated burns is three dimen- 
sional and uses the optimum steering law. 

The analysis presented herein assumes that all nodal correction is done during the 
outbound burn and the retroburn. An alternative to this approach would be to introduce 
an additional nodal correction burn, or some nodal correction could be done during the 
impulsive apogee and perigee burns. These approaches were not analyzed, and it is not 
clear if they offer any theoretical improvement over the approach used in this report. 
Based on the results obtained, the added nodal correction burn, even if theoretically 
optimum, may not be justified since the payload loss due to nodal correction is quite 
small as is shown in the RESULTS AND DISCUSSION section. Although small payload 
gains might be realized by the addition of a nodal correction burn, the losses associated 
with the additional startup and shutdown of a real engine could substantially reduce the 
expected payload gain. 

In addition to these constraints the following assumptions were made in deriving the 
equations of motion governing the problem: 

(1) A spherical nonrotating Earth model is assumed. This implies that there is no 
nodal or apsidal precession of the Tug orbits and that the Tug will correct for nodal pre- 
cession of the SS orbit only. This gives a somewhat conservative vehicle performance 
since, if the Tug orbits were allowed to process, the total nodal correction required 
would be slightly less and consequently the payload capability would be greater. For 
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the sample cases analyzed, nodal precession of the Tug orbits is less than 0. 2° com- 
pared with a total SS-orbit nodal precession of approximately 8°. 

(2) The right ascension of the outgoing asymptote is not constrained in the trajectory 
optimization. Instead an optimal "pseudo right ascension" is obtained for each declina- 
tion of the outgoing asymptote. Pseudo right ascension is defined as the longitude of the 
projection of the outgoing asymptote in the equatorial plane measured from the ascending 
node of the initial SS orbit. The true right ascension, then, is the sum of the pseudo 
right ascension and the longitude of the SS-orbit ascending node as measured from the 
vernal equinox. The longitude of the SS-orbit ascending node is a function of the SS 
launch time and launch azimuth. For a fixed SS launch azimuth and a 1-day launch win- 
dow the longitude of the SS-orbit ascending node sweeps through 360°, and any desired 
right ascension can be obtained simply by selecting the proper SS launch time. 

(3) The reusable Tug and the expendable kickstage are assumed to have a constant 
engine exhaust velocity (specific impulse) and thrust, although the analysis is formulated 
to accommodate a variable thrust (and flow rate). 

(4) The analysis is done for a two-body problem, and all perturbations due to the 
Sun, Moon, and planets are neglected. 

(5) The Tug departs from and returns to a circular SS orbit. The longitude of the 
line of nodes of the initial SS orbit is fixed, and the Tug is free to depart at an optimum 
point along this orbit. The assumption of a circular SS orbit greatly simplifies the 
return-orbit geometry and rendezvous requirements. For elliptic SS orbits the line of 
apsides of the final Tug orbit must coincide with that of the precessed SS orbit, and addi- 
tional constraints would have to be imposed on the problem. This apsidal constraint is 
discussed in more detail in the following section. 


Nodal and Apsidal Precession 


Nodal and apsidal precession were analyzed and analytic equations were derived in 
reference 3. For convenience they are reproduced here. Symbols and notation used in 
the following equations and the analysis are defined in appendix A. 


Ail 


jVg 




2 3/2 

(1 - e^) cos i 


( 1 ) 


and 
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A(t> = 


7/2 


( 2 ) 


2p|/2 \p 


2 2 
(1 - e^) (4 - 5 sin^i) 


Nodal precession AS2 and apsidal motion Aw are functions of the inclination, eccen- 
tricity, and semilatus rectum of the particular orbit. They also depend on the total 
time Tj^ spent in orbit (total trip time). In equations (1) and (2), J, G, and are 
the second harmonic in the gravitational potential function (oblateness parameter), the 
gravitational constant, and the equatorial radius of the Earth, respectively. 

Equation (1) is used in the analysis to compute SS- orbit nodal precession. Note that 
the line of nodes moves at a constant rate in the retrograde direction for the given SS 
orbit. 

Apsidal motion, given by equation (2), becomes an important consideration if ellip- 
tic SS orbits are used. For circular orbits the apsidal precession translates into an 
effective reduction in the period of the SS orbit. However, this can easily be taken care 
of with minor payload loss or gain by varying the total trip time a few minutes. Al- 
though theoretically proper phasing of Tug and SS can be achieved by varying total trip 
time, the operational aspects of this approach on Tug guidance hardware and software 
must be carefully evaluated. These requirements are outside the scope of this report 
and are not considered further. 


Equations of Motion Governing the Problem 

Based on the assumptions discussed in the previous sections equations of motion 
governing the problem are presented here. The variables and other notation used in 
this report are defined in appendix A. Equations governing the flight of a rocket in an 
inverse-square gravitational field are given by 


G V <3(0 - 

- — r + ® f 

J.2 m 

(3) 

? = V 

(4) 

I 

II 

• a 

(5) 

k = 1 

(6) 

f • f - 1 = 0 

(7) 
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In these equations is the rocket exhaust velocity, |3 the mass flow rate as a function 
of time, and f the unit thrust direction. The superscripts “ and ' represent vector 
and unit vector quantities, respectively, and • denotes the time derivative of the par- 
ticular variable. State variables r, V, m, and ^ are radius, velocity, mass, and a 
dummy variable associated with mass flow rate, respectively. By using a Hamiltonian 
formulation of the Maximum Principle, adjoint equations associated with the state vari- 
ables can be obtained from the Hamiltonian given by 


H = X • V + M • r + om + T? + y(f • f - 1) (8) 

The terms X, II, a, and t in equation (8) are the adjoint variables associated with the 
state variables V, r, m, and respectively. The thrust magnitude constraint given 
by equation (7) is adjoined to H by y. 

From equation (8) the time derivatives of the adjoint variables on each trajectory 
phase are given by 


where 


X - -II 


M = [x - 3(X • f)f 1 

3 

r 


_ 


a = 


(A • f) 


m 


a? 


V 

K = (X • f ) - a 
m 

The thrust direction that minimizes equation (8) is given by 

f =X 


(9) 

( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 


and 
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y = - 


(15) 


2m 

To solve the multipoint boundary value problem determined by some initial, final, and 
intermediate boimdary conditions, state equations (3) to (5) and adjoint equations (9) to 
(12) must be numerically integrated. Optimum thrust direction is along the unit X vec- 
tor given by equation (14). 


Boundary Conditions and Constraints 

In the following paragraphs, problem constraints and the mission profile are dis- 
cussed in more detail. Sketch (b) illustrates the simplified trajectory profile with the 


/.^Outgoing 
/ asymptote 

/ 

/ 

/ 



Tug and expendable Mckstage departing from a given initial orbit. Outbound Tug burn 
is initiated at time a^ = 0 at an optimum argument of latitude in the initial circular 
inclined orbit. Thrust is terminated at time bj^, followed by a fixed coast phase to time 
b 2 . During the coast phase, the payload and kickstage (payload attached to the kickstage) 
are separated from the Tug. The Tug is turned around, and a safe separation distance 
is established for both kickstage firing and Tug retroburn firing. The Tug retroburn 
and the fixed kickstage burn start at times a^ and a^ (a^ = a^ = b 2 ), respectively. At 
kickstage burnout b^ the desired energy and declination must be achieved. Tug retro- 
burn is terminated upon reaching an orbit that will achieve the specified total trip time 
when followed by apogee and perigee delta velocity maneuvers. With these time conven- 
tions, the constraints governing the problem can be written mathematically as follows: 

qj = r(aj) - r j = 0 (16. 1) 
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Q2 = V(ai) - Vi = 0 

(16.2) 

qg =r(aj) • V(aj) = 0 

(16.3) 

= ii(aj) • 2g - cos ij = 0 

(16.4) 

= -h(aj) • ?2 ■ ^1 ij = 0 

(16.5) 

qg =m(aj) - = 0 

(16.6) 

«l7 = C(aj) - aj = 0 

(16.7) 

qg=ai =0 

(16.8) 

qg = r(a 2 ) - r(bj) =0 

(16.9) 

qiO = V(a2)- V(b^) = 0 

(16. 10) 

q^j = m(a2) - ni(bj^) = 0 

(16.11) 

q^2=?(a2)-?(bi)=0 

(16.12) 

qjg =r(ag) - r(b 2 ) = 0 

(16.13) 

qi4=V(ag)- V(b2) = 0 

(16.14) 

q^g = m(ag) + mCa^) - m(b 2 ) = 0 

(16.15) 

q^g = ^(ag) “ ^(^ 2 ) ~ ^ 

(16.16) 

qi7 =r(a4) - r(b2) =0 

(16.17) 

qi8=V(a4)- V(b2)=0 

(16.18) 


qi9 = C(a4) - ^(bg) = 0 

(16.19) 

^20 = ag + At - bg = 0 

(16.20) 

q2i =f(bs) • Hq = 0 

(16.21) 

^22 = V(bg) • = 0 

(16.22) 

1 

q 23 = m(bg) - m^e ’ = 0 

(16.23) 

^124 = ’’3 - ^1 + = 0 

(16.24) 

q25=E[r(b4),V(b4)]-Ej5 = 0 

(16.25) 

q26 = <;p[r(b4),V(b4)]- <pj^ = 0 

(16.26) 

qg7 = m(a4) - m(b4) - mpj^ = 0 

(16.27) 

^128 = a2 - = 0 

(16.28) 

‘^29 = ^3 " ^2 ^ ® 

(16.29) 

^30 = a4 - bg = 0 

(16.30) 


Equations (16. 1) to (16. 8) describe the given initial conditions at the beginning of 
outbound Tug burn. The only free orbital element is the argument of latitude where the 
burn starts, which will be used as an initial condition in solving the two-point boundary 
value problem. Equations (16. 9) to (16. 20) are constraints expressing the continuity of 
state variables at times a. 2 > 2 .^. In equation (16. 15), m(a^) is the combined 

kickstage and payload mass and, in equation (16. 20), At is a fixed coast time. The 
Tug is constrained to return to a desired final orbit that coincides with the precessed SS 
orbit at rendezvous. This constraint is expressed by equations (16. 21) and (16. 22). 

Tug retroburn is terminated when the total-trip-time constraint (eq. (16. 24)) is satisfied, 
and at this point the Tug must have sufficient propellant remaining for the ideal impul- 
sive apogee and perigee maneuvers required to return the Tug to the SS (eq. (16. 23)). 

The kickstage after a fixed burn (eq. (16. 27)) must satisfy a desired energy (eq. (16. 25)) 
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and declination of the outgoing asymptote (eq. (16. 26)). Equations (16. 28) to (16. 30) 
establish continuity in time at staging points along the trajectory. By using the theory 
of branched trajectory optimization (refs. 4 and 5) a set of variational boundary condi- 
tions are derived in the following paragraphs from the constraints previously described 
and an augmented Hamiltonian given by 


H* =-m(b4) + 2]Hi*|^€jqj 
i=l i=l 

where m(b^) is the mass at kickstage burnout that will be maximized and are the 
Hamiltonians during outbound burn, Tug/kickstage separation coast time, Tug retroburn, 
and kickstage burn. Constraints qj. are adjoined to H* by usir^ the arbitrary multi- 
pliers e. . Auxiliary variational conditions associated with the multipoint boundary value 
problem can be derived by using this augmented Hamiltonian. 


Auxiliary Variational Equations 

The constraints described in the previous section and the augmented Hamiltonian 
are used to satisfy the following auxiliary variational equations in order to maximize 
payload : 


X(a^) = + e 3 r(a^) + V^q^ + V^q^ 

(17.1) 

ju(a^) = e^f(ap + e3V(a^) + ^q^ + V-q^ 

(17.2) 

^(^ 2 ) - ^(bj^) 

(17.3) 

Mag) = Mb^) 

(17.4) 

a(a2) = cr(b2) 

(17.5) 

T(a2) = r(b^) 

(17.6) 

^(ag) +X(a4) = X(b2) 

(17.7) 

Mag) + Ma4) = ju(b2) 

(17.8) 
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II 

(17.9) 

T(a^) + r(a^) = T(b2) 

(17.10) 

^(a^) = a(b2) + a(b^) - 1 

(17.11) 

. m(b^) 

X(b ) + — ^ h^ + a(b„) ^ AV.p + H(b„) V„T = 0 

^ V(b3) ° Vg^^ ^ ^ ^ 

(17.12) 

^ m(bJ 

M(bJ + h + a (b„) ^ V- AV^ + H(b ) V-T = 0 

r(b 3 ) ° Vg 1 ^ ^ 3 r 

(17.13) 

o 

II 

(17.14) 

X(b4) + €25 ^V^ ^26 ^V*^ " ° 

(17.15) 

M(b^) + €25 ^r® ^26 ^r^ " ° 

(17.16) 

T(b4) = 0 

(17.17) 

H(b4) = 0 

(17.18) 

H(a2) + H(a3) + H(a4) - H(b^) - H(b2) = 0 

(17.19) 


In these equations, are arbitrary constants that must be evaluated to obtain the 
desired variational final conditions. Variational boundary conditions at the beginning of 
outbound Tug burn are given by equations (17.1) and (17.2). There are six of these 
equations and five arbitrary constants. Therefore, there is one variational condition 
that can be used as a final condition; it will be derived later. Equations (17.3) to (17.6) 
give the continuity of adjoint variables at the end of the outbound Tug bxirn; equations 
(17.7) to (17.11) give the conditions at the end of the Tug/kickstage separation coast 
phase. Variational conditions at retroburn completion are given by equations (17.12) to 
(17.14), and variational conditions at kickstage final bvu:nout by equations (17.15) to 
(17.18). The Hamiltonian evaluated at various times along the trajectory must satisfy 
equation (17.19) . 
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Gradients of return time T, total delta velocity AV>p, and declination cp with 
respect to r and V have been derived in reference 1 and will not be reproduced 
herein. 

The variational auxiliary conditions given by equations (17) must be satisfied by the 
optimal trajectory in order to maximize payload. The arbitrary constants in these 
equations and the available variational boundary conditions are evaluated in the following 
sections. Conditions at the end of outbound Tug burn and at the start of Tug/kickstage 
separation coast phase can all be satisfied directly (eqs. (17.3) to (17.6)) by setting the 
corresponding variables at ag equal to those at b^ and will not be discussed in more 
detail. 

As has been mentioned, some of the auxiliary variational equations contain arbi- 
trary constants that must be evaluated or eliminated from the equations to obtain the re- 
quired variational boundary conditions. For example, there are six equations given by 

(17.1) and (17.2) with five arbitrary constants. The arbitrary constants could be eval- 
uated from any combination of five of these equations. Once the constants are known, 
the resulting independent equation can be used as a free boundary condition. The form 
of this boundary condition will depend on the choice of eqxiations used in computing the 
arbitrary constants. An alternative to computing these arbitrary constants is to find 

a function (a combination of two three-dimensional vector equations, such as (17.1) and 

(17.2) , and the state conditions) that is not a function of these arbitrary constants and to 
use this as the variational boundary condition to be satisfied. This, of course, assumes 
an a priori knowledge of the function to be used. Such a function has been found and will 
be used to evaluate some of the variational boundary conditions. When the form of this 
function is not known, the arbitrary constants will be evaluated. And the variational 
boundary conditions will be computed in a convenient (not necessarily orthogonal) coor- 
dinate system to assure that the variational conditions are independent. 


Variational Condition at Departure from Initial Orbit 

As was mentioned previously, there are six independent variational equations (eqs. 
(17.1) and (17.2)) that must be satisfied at departure from the given initial SS orbit. 
This set of equations has five arbitrary constants, leaving one variational condition that 
is available as a boundary condition. The free variational condition can be derived as 
follows: Define 


C(t) = X(t) X V(t) + /i(t) X r(t) (18) 

where x represents the vector product of two vectors. From equations (17. 1) and 

(17 . 2) , the C vector at time a^ is given by 
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Substituting the gradients of q^ and q^ with respect to V(aj^) and r(aj^) into C(aj^) 
gives 


C(a.j^) — (^4^2 “ ^5^2^ ^ ( 19 ) 

Equation (19) implies that C(a^) is perpendicular to li(aj^); consequently, their scalar 
product is zero, which is the required free variational condition 

C(aj) . h(a^) = 0 (20) 

It is shown in reference 6, that the free initial adjoint variables X(aj^) and can be 

replaced by a set of physically more meaningful parameters; namely, pitch attitude if/, 
pitch rate 1^, yaw attitude 6, yaw rate 5, and magnitude of X and X. The derivatives 
of X can be computed from equation (20) as shown by the same reference . 


Variational Conditions at Kickstage Burnout 

In this section, variational conditions are derived at kickstage burnout. The arbi- 
trary constants egs ^26 equations (17.15) and (17.16) can be evaluated as fol- 
lows: Compute C(b^) from equation (18), giving 

C(b4> = -eggy (21.1) 


where 


y = X V(b4) + Vpp X r(b^) (21 . 2) 

The magnitude of vector y is equal to 1 (as can be shown by using eq. (B5)); and from 
equation (21 . 1) , €20 is given by 


^26 = • y (22) 

Vector y is perpendicular to = 0,0,1) as shown in appendix B, and a free varia- 

tional condition may be obtained from equation (21 . 1) by 
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CCb^) • ^3 = 0 


(23) 


The constant egs chosen to satisfy the magnitude of equation (17.15); that is, 


"25 = ± 


|Mb4) + €26 
VCb^) 


(24) 


The choice of sign is determined as follows: Suppose the declination is arbitrary (then 
€26 = 0)' With this assumption after substituting €25 ^26 (l*^* 

that equation becomes 


X(bJ _ 

Mb.) =T_l.V(b4) 
^ V(bJ ^ 


But the thrust direction is along X and, at kickstage burnout, energy must be added es- 
sentially along the velocity vector in order for a desired energy to be reached. This 
indicates that the minus sign should be used in equation (24). 

Another free variational condition is obtained by substituting equations (17.15) and 
(17.16) into equation (8) 


H (» 4 ) - -‘26 


— . r(b^) + Vpp . V(b^) 


+ ^(^4) 


■ ■ - <r(b4) 


m(b.) 


+ T(b4) 


However, 


- ^Y<P • r(b4) + . V(b4) = ^ 

r^(b4) 


along a coast phase, but declination along a coast phase is constant and d<^/dt = 0, 
giving 


H(b4) = ^(b4) 


Ve.2^(^4) 


m(bA 


- <^(b4) 


+ T(b4) 


Substituting equations (17.17) and (17.18) into this equation gives 
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|3(b4) 


Ve.2^(b4) 

mCb^) 


- crCb^) 


= 0 


In this equation, /Slb^) + 0 since, in order to reach the desired energy, thrust must be 
on and therefore. 


crib.) = ± ~ 

^ m(bj 


Substituting this equation and equations (17.17) and (17.18) into equation (8) at time b^ 
gives the desired variational boundary condition 


. __G_X(b 4 ) • ^(^4^ + ^(^ 4 ) • V(b 4 ) = 0 (25) 

r^b4) 


Two additional variational final conditions can be derived by taking the scalar product of 
equations (17.15) and (17.16) with V(b^) and r(b 4 ), respectively, giving 


[x(b4) + €25V(^4^ + ^26 ' ^^^4^ = ° 


(26.1) 


M(b4) + €35 


^ r(b 4 ) + £26 


r^(b4) 


■ ?(b4) = 0 


(26.2) 


Thus, there are four variational conditions at b 4 that can be used to solve the two- 
point boundary value problem given by equations (23), (25), (26.1), and (26.2). 


Variational Conditions at Tug Retroburnout 

The six equations at retroburnout given by equations (17.12) and (17.13) have three 
unknowns, the two arbitrary constants f21>^22 °'(b 3 ). To determine the three 

variational conditions implicit in the six equations, these unknowns must be determined. 
The arbitrary constants Cgi ^22 are evaluated from C(bg), given by 

C(b3> = \ X X €2if(bj)] (27) 
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This equation was derived from the results of reference 7, where it is shown that 
VyX X V + VjX X r = 0 if X = x[r, V,r • v]. Since total delta velocity AV,p and return 
time T are of this form, the corresponding vector products are zero. Equations (27) 
implies that C(b 2 ) is perpendicular to h^ and that their scalar product is zero. One 
free variational condition is given by 


C(b 3 ) • ho = 0 (28) 

The constants and 622 are evaluated by taking the scalar product of equation (27) 
with v{h^) and ^(bg), giving 


«b3> • l>(bj) ^ 


(29.1) 


X(b,) . h(b„) ^ 

'22= ^V(bg) 

22 . h(b3) 3 


(29 . 2) 


The adjoint multiplier cr(b 3 ) is evaluated by substituting equations (17.12) to (17.14) 
into equation ( 8 ): 


H(b„) 


1 - 


r^bj) 


VyT . r(b3) + • V(b3) 




m(b„) 


e,l 


r^(b,) 


Vy AV^ • r(b3) + V-AV^ • V(b3) 


+ ^(bg) 


m(bJ 


CT(bJ 


(30) 


However, 


G 

r2(b3) 


V^T . r(b 3 ) + V-T • V(b 3 ) = ^ 
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and 




d AV, 

" V(b3) 

^ dt 


These equations are valid on a coast phase. Ho-wever, the change in return time is equal 
to the change in time along the coast arc but opposite in sign (dT/dt = -1), and the delta 
apogee and perigee velocities required to return the Tug remain constant along a coast 
arc (d AV,p/dt = 0) . Substitutii^ these results into equation (30) gives 


^(bg) 




m(b,) 


- 


= 0 


The mass flow rate ^(bg) cannot be zero since thrusting must continue until a desired 
total trip time is reached at time bg. Therefore 


a(bg) = 


Ve.l"(b3) 

m(b„) 


(31) 


Two variational conditions can now be derived by taking the scalar product of equa- 
tion (17.12) with AV,p and the scalar product of equation (17.13) with and 

substituting u(bg), €2^, and €22 


X(bg) • h(bg) ^ 

X(bg) + X(bg) AVt + H(bg) V^T 

ho • h(bg) 


V^AV,p = 0 (32.1) 


M(bg) • h(bg) ^ 

/Kbg) . X(bg) V- AVt - -J-^h^ H- H(bg) V-T 

% ■ hCbg) 


V- AV^ = 0 


(32.2) 


where 


H(bj) = - _£_x(b3) • rCbj) + ^(bj) ■ VCbj) 

r=(b3) 

It is important that these derived variational conditions be independent, that is, that 
there exist no trivial solutions or that conditions previously derived do not already 
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satisfy these equations. It is easy to see that conditions (32. 1) and (32, 2) are inde- 
pendent by looking at the terms ^(bg) | AV,j, 1 ^ and A(bg) | V- AV,p | ^. These terms 
cannot be zero since the total delta velocity to return the Tug cannot be zero. The 
gradients cannot be zero since any change in radius and velocity will change AV,j.. 
Also, the magnitude of A(bj) is a function of these gradients and hence cannot be zero. 
Since these terms cannot be zero, a trivial solution cannot exist and the conditions are 
independent. 

The three free variational boundary conditions derived herein are given by equa- 
tions (28), (32.1), and (32.2) 


Variational Conditions at Tug/Kickstage Separation 

Equations (17.7), (17.18), and (17.10) express constraints that the adjoint variables 
must satisfy on an optimum trajectory . Adjoint variables X and JI are known at time 
bg (they are integrated from time a^^) . If these variables were known at time ag, the 
preceding equations could be used to compute them at a^. At time ag, X and JI are 
arbitrary; therefore, ^//(ag), i^(ag), 6(ag), 6(ag), X(ag), and x(ag) can be used as initial 
conditions to compute X(ag) and ju(ag) . The derivative X(ag) is computed from equa- 
tion (28) by using the procedure given in reference 6. The adjoint variable X(ag) can 
also be computed as follows: From equations (17.7) and (17.8) the relation between the 
C vectors at kickstage separation is given by 

C(ag) + C{n^) = C(b 2 ) (33) 

The vector C is a constant of motion on every subarc (C = constant) . This can be veri- 
fied by taking the time derivative of C. After substituting for X, /u, r, and V, the 
time derivative C becomes identically zero, which implies that C is a constant. 

Since C is a constant, C(a^) = C(b^). Therefore, taking the scalar product of this 
equation with Zg and substituting equation (23) give C(a^) • ig = 0, and X(ag) can be 
computed from equation (33) as 


X(ag) = 


C(b2> • 


C(aj) 


x(ao) 


(34) 


where 
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C(aJ 


^(ag) 


r^(a<}) ^ V(ag) “ ^(ag) x r(ag)"] * h x -i 

= X(ag) X V(ag) |"C/„ \ \1 . u ^ ^ ^^^3^ 


|X(a 3 ) X r(ag)] • 


the vectors Ma^) and X(ag) are evaluated as sho"wn in reference 6. Therefore, equa- 
tion (34) can be used to compute the magnitude of X at ag. 

Another variational final condition can be derived from equation (17.19). 

H(a 2 ) + H(ag) + H(a 4 ) - H(bj) - H(b 2 ) = 0 
Since H is constant on every branch, 


H(a2) = H(b2) 


And from equations (17.3) and (17.4) 


^ - • r(b 2 ) + ju(b 2 ) ‘ V(b 2 ) = - X(bj) • r(b^) + /I(b^) • V(b^) 


r^(bJ 


r^(b,) 


Also from equation (17.18) and the constancy of H 


H(a4) = H(b4) = 0 

Substituting these results gives 

H(ag) = H(bi) (35) 

From equation (11) it is obvious that a is constant' on a coast phase, and from 
equations (17 . 5) and (17 . 9) 


a(ag) = CT(bj^) 

Without loss of generality, assume /3(bj) = ^(Sg). Substituting these results and equa 
Won (17.10) into equation (35) gives 


— ^-(a^) • r(a^) + V(a 4 ) • M(a^) + T(a^) + ^)3(b^) 

r^(a^) 


A(bj) 

m(b^) 


m(ag) 


(36) 
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If the kickstage flow rate is constant, the adjoint variable t also remains constant 
during the burn, and from equation (17,17) T(a^) =0, If the kickstage flow rate is vari 
able, r(a^) is computed from 


K M dt 


(37) 


Initial and Final Conditions 

Some of the initial conditions have been discussed in previous sections. In this sec- 
tion a complete set of initial and associated final conditions are given: 


Initial conditions 

Final conditions 



?(bg) • ho = 0 

(38.1) 


o 

II 

o 

<> 

(38.2) 

6(a^) 

Equation (26. 1) 

(38.3) 

6(a^) 

Equation (26.2) 

(38.4) 

X(aj) 

Equation (17.11) 

(38.5) 

Maj) 

Equation (20) 

(38.6) 


m(a^) - m(b^) - m^^ = 0 

(38.7) 

V^Cag) 

- <Pd = 0 

(38.8) 

^(ag) 

Equation (25) 

(38.9) 

6(ag) 

Equation (32. 1) 

(38.10) 

Mag) 

Equation (32. 2) 

(38.11) 

X(ag) 

Equation (23) 

(38.12) 
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Equation (28) 


(38.13) 


Mag) 


b 


3 


b 


4 


m(ag) 

u(a^) 



+ T - Tp = 0 


E(b^) -Ej5 = 0 

1 

m(bg) - nij^e ’ = 0 

Equation (36) 


(38.14) 

(38.15) 

(38.16) 

(38.17) 


In solving the two-point boundary value problem, the 17 initial conditions must be 
guessed at and a numerical iteration process performed vmtil the given final conditions 
are satisfied. The iteration size can be reduced to an 11-by-ll iteration as follows: 
Terminate the Tug retroburn and the kickstage burn when equations (38.15) and (38.14) 
are satisfied, respectively. Furthermore, equations (38.6), (38.12), and (38.13) can 
be used to compute A.(aj), A(ag), and A(ag), respectively; and final condition (38.5) can 
be satisfied trivially by scaling the problem with A(a^) . 

The numerical results presented were obtained by using a Newton-Raphson iteration 
technique. Each initial condition was perturbed individually and the changes in final con- 
ditions observed in order to evaluate the partial derivatives of the final conditions with 
respect to the given initial conditions. The perturbation sizes of the initial conditions 
were determined empirically. Through trial and error it was determined that 0.1 per- 
cent perturbation of the initial conditions gave sufficiently accurate partial derivatives. 
Using these partials showed convergence for the sample cases to be well behaved and 
uniform. Some convergence difficulties were encountered for declinations whose abso- 
lute value exceeded initial SS orbital inclination. For these cases, payload decreases 
very rapidly with increasing declination of the outgoing asymptote, and the finite - 
difference partial derivatives become very sensitive to initial conditions. The numer- 
ical iteration process was continued until the expected percentage of payload improve- 
ment and the normalized final conditions were within a given tolerance. 


RESULTS AND DISCUSSION 

The method and equations developed in the earlier sections were applied to a given 
vehicle configuration and a range of planetary target conditions . The results of this 
analysis are presented in this section. 

Nominal mission constraints and planetary injection conditions used in the sample 
case analysis are given in table I. As is discussed in the INTRODUCTION, a single - 
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stage Tug in the reusable mode can only be used for relatively low-energy missions. 

To extend the energy range of this reusable Tug to the region covering comet and outer 
planet missions an expendable kickstage was introduced. The reusable Tug delivers the 
kickstage and payload to some energy beyond Earth escape, and the kickstage provides 
the additional energy required for the mission. The energy range considered for the 
baseline vehicle configuration is from a vis-viva energy of 35 to 120 km /sec . The low 
end of the energy range is determined by the performance of the single-stage reusable 
Tug. At the high end the payload capability of the present configuration approaches 
zero. Departure -orbit altitude and inclination were chosen to be 185 kilometers and 
28.5°, respectively. This represents a due-east SS launch from the Eastern Test 
Range, Based on previous experience with single-stage reusable Tug trajectories 
(ref. 1), the total trip time was chosen to be 1 day and Tv^/kickstage separation coast 
time was selected as 10 minutes. Althov^h the 1-day trip time is nearly optimum for 
performance, the selection of 10 -minute separation time was arbitrary, and significant 
payload gain can be realized if this time can be shortened. 

In the analysis it is assumed that the Tug, kickstage, and payload are launched with 
a single SS. This assumption implies that the weight of the selected configuration is 
fixed (28 622 kg) and consistent with SS payload capability. To keep the initial weight 
fixed, the Tug propellant must be varied to accommodate changes in payload capability. 

Baseline Ti^ characteristics are given in table II. Engine performance corresponds 
to that of an uprated RL-10 with a specific impulse of 460 seconds and an engine thrust 
of 66 723 newtons. The RL-10 engine is currently in use on the Centaur stage. 

The characteristics of the space -storable kickstage were taken from reference 8 
and are given in table III. The kickstage is pressure fed and has a specific impulse of 
378 seconds. Because of its low thrust level (13 345 N) and the spatial location of this 
burn, the kickstage can be expected to have large gravitational losses. 

As a first step in this study the range of useful payload was determined for this 
vehicle configuration. Maximum payload capability is determined as a function of mis- 
sion energy for a 0° declination. Selection of this declination is arbitrary, and it is 
intended to establish the useful payload range of this vehicle configuration as a function 
of mission energy. These results are presented in figure 1. For comparison, ideal 
impulsive payload capability without nodal correction is also given. Ideal impulsive 
performance was computed by assuming an infinite thrust for the Tug and kickstage and 
that all velocity change maneuvers would take place instantaneously at the same low al- 
titude. Ideal, impulsive, mission delta velocities were computed based on the given 
mission energy. All velocity changes were assumed to be planar, and they do not ac- 
count for plane change requirements due to nodal precision or changes in inclination to 
meet mission declination requirements. Note the large payload losses. At an energy 
of 100 km /sec the ideal payload is 1200 kilograms and the loss due to gravitational 
and thrust -to -weight effects, nodal precession, etc., is 520 kilograms. Obviously, 
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ideal impulsive performance calculations for this configuration and mission are inade- 
quate, and the trajectories must be integrated to include actual trajectory losses. To 
further analyze the behaAaor of these trajectories two energies were selected; an en- 
ergy of 50 km /sec , which corresponds to a high-energy inner planet or comet mis- 
Sion, and an energy of 100 km /sec , which corresponds to an outer planet mission. 

To explain the results that follow, it is helpful to recall some facts from orbital 
geometry and previous analyses of planetary missions (ref. 1). If the SS-orbital nodal 
precession is assumed to be zero and a spherical Earth model is used, for a fixed SS 
orbit the same payload can be delivered, using planar trajectories, to all declinations 
of the outgoing asymptote (DLA) whose absolute value is less than the initial SS orbital 
inclination without using yawed trajectories. Furthermore, there are two distinct solu- 
tions for each DLA whose absolute value is less than the SS orbital inclination. As the 
DLA approaches the SS orbital inclination the two solutions approach one another, and 
at a DLA equal to the SS orbital inclination they degenerate into a single solution. Those 
DLA's whose absolute value is greater than the initial SS orbital inclination can be 
reached only if the Tug orbital inclination is increased to or beyond the absolute value 
of the desired DLA. To change orbital inclination at low -Earth -orbit altitudes is very 
costly, therefore, maximum performance for these DLA's is obtained when the orbital 
inclination is made equal to the absolute value of the desired DLA. 

Another important point to recall is that, for a fixed out -of -plane acceleration, 
maximum change in the line of nodes is obtained when the out-of -plane acceleration is 
applied at the antinode and maximum change in orbital inclination is obtained when the 
yaw maneuver is done at the line of nodes (at the equator). An out-of -plane accelera- 
tion, which tends to increase inclination near the line of nodes, will move the line of 
nodes retrograde when applied before the nodal crossing and posigrade when applied 
after the nodal crossing. Obviously, since the line of nodes of the Tug orbits must be 
moved in a retrograde direction, inclination changes that move the node in this direction 
are desirable. An inclination change that moves the node posigrade is undesirable be- 
cause it increases the total nodal correction that must be made later . 

Based on equation (1) and nominal trajectory and mission parameters, the SS 
orbital precession is approximately 8° in a retrograde direction. The Tug will have to 
correct for this nodal precession during outbound burn and retroburn. Since the plane 
change maneuvers during Tug outbound burn and retroburn are small, the basic trajec- 
tory profiles are not expected to be drastically different from the planar trajectory pro- 
file without nodal correction. And the general conclusions drawn from the planar case 
are expected to apply to the more complex nonplanar trajectories. 

Existence of the two distinct solutions for DLA's whose absolute value is less than 
the SS orbital inclination was shown numerically for the present problem. This was 
done by taking a converged case for a given DLA, changing the signs on initial yaw alti- 
tudes and rates, and adding 180° to the Tug departure argument of latitude. With these 
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initial conditions the trajectory was integrated and without any iteration gave a con- 
verged case for a DLA of the same magnitude but opposite sign. Both solutions will be 
presented for the two mission energies, although for simplicity the discussion will be 
restricted to one of the solutions. To eliminate confusion, results obtained from the 
second solution are shown by a dotted line in figures 2 to 13. Discussion of the results 
will be restricted to the first solution (solid line), but similar reasoning can be applied 
to the second solution. 

Data presented in figures 2 to 7 as a function of the declination of the outgoing 
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asymptote are for the 50-km /sec energy case. Figure 2 gives the optimum pseudo 
right ascension. As was discussed earlier (in the section Flight Profile and Simplifying 
Assumptions), specific requirements on the right ascension of the outgoing asymptote 
can be satisfied for a given optimal pseudo right ascension by selecting the proper SS 
launch time . By subtracting the value of pseudo right ascension for the desired DLA 
from the required right ascension, the longitude of the ascending node for the SS orbit 
may be determined. Since the SS orbital longitude of the ascending node sweeps through 
360° for a fixed SS launch azimuth and a 1-day launch window, any desired right ascen- 
sion of the outgoing asymptote may be obtained by selecting the proper SS launch time . 

The point of Tug departure from the initial orbit was optimized and the results are 
presented as a function of DLA in figure 3. Departure argument of latitude is defined 
as the angle measured in the initial SS orbital plane from the ascending node of this 
orbit to the point where the outbound burn starts. Note the similarities between this 
and the previous figure. To explain the similar behavior of these figures from orbital 
geometry, recall that the shape of the escape hyperbola depends primarily on perigee 
radius (which changes little with DLA) and mission energy (which is constant). There- 
fore, once the trajectory is optimized for a given DLA, the shape of the escape hyper- 
bola does not change significantly as DLA is varied. As DLA is changed, the plane of 
the escape hyperbola is essentially rotated about the angular momentum vector of the 
initial SS orbit until the desired DLA is satisfied. Therefore, as departure argument of 
latitude decreases so must pseudo right ascension. 

Orbital inclination at Tug/kickstage separation and nodal correction during outbound 
Tug burn are presented in figures 4 and 5, respectively. To explain the trajectory 
characteristics obtained, recall the earlier discussion on the effect of out -of -plane 
maneuvers on the orbital inclination and the line of nodes. For negative DLA's whose 
absolute value is less than the SS orbital inclination the outbound burn starts past the 
line of nodes (fig. 3, solid line) and continues through a burn arc of approximately 100°. 
Thus, a major protion of the outbound burn occurs near the antinode, which is an opti- 
mal region for changing the longitude of the line of nodes. Almost all the required 
nodal change (8°) is done during this burn, as seen in figure 5. For negative DLA’s 
whose absolute value is larger than the SS orbital inclination the orbital inclination must 
be increased to reach the desired DLA, as seen in figure 4. Orbital inclination at the 
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completion of the outbound Tug biurn is less than the absolute value of the corresponding 
DLA. Since final spacecraft orbital inclination must be equal to the absolute value of 
the DLA, as was discussed earlier, the remaining inclination change is done during the 
kickstage burn. Note that nodal change during the outbound burn decreases as the ab- 
solute value of the DLA is increased beyond the SS orbital inclination. This is caused 
by the large inclination change being made during the outbound burn. Since the Tug must 
return to the SS orbit, the same inclination change must be made during retroburn, and 
it becomes more efficient to perform a larger nodal correction during this burn. 

As the value of the DLA is increased from negative values to zero, the departure 
argument of latitude moves closer to the line of nodes and it becomes more difficult to 
make large nodal corrections during the outbound burn (fig. 5) . Also for these DLA's 
the burn takes place past the line of nodes and orbital inclination must be decreased to 
move the node retrograde, as seen in figure 4. For DLA's between 10° and 20° the 
outbound burn traverses the line of nodes and nodal change for these cases is a minimum, 
as expected. As the DLA is further increased the outbound burn begins well before the 
nodal crossing and the nodal correction increases. Also orbital inclination increases 
to move the node retrograde. For DLA's larger than the SS orbital inclination a large 
inclination change is being made during the retroburn, and again it becomes more effi- 
cient to make larger nodal corrections during this burn. Consequently, nodal correc- 
tion during outbound burn (fig. 5) decreases. Note the single solutions for high negative 
and positive DLA's. The reason these curves do not leave the closed curve exactly at 
DLA's equal to the SS orbital inclination is that orbital inclination must be changed dur- 
ing the outbound burn to correct for nodal precession. 

Payload as a function of DLA is given on figure 6. Payload is maximum for DLA's 
of approximately ±12°, and it varies less than 1 percent for DLA's between -28. 5° and 
+28. 5°. Payload outside this range decreases rapidly because of the large inclination 
changes being made. Note that payloads for the two solutions are very close althov^h 
the trajectory profiles are different. This is a result of the trade-off between losses 
during outbound burn and retroburn. 

Velocity losses during the three main trajectory phases (outbound Tug burn, return 
leg (inbound), and kickstage burn) are given in figure 7. Velocity loss is defined as the 
velocity increment supplied by the vehicle (computed from the ideal rocket equation and 
the amount of propellant used by the vehicle) above the ideal impulsive velocity (as de- 
fined earlier) required to perform the mission. By this definition, the velocity loss 
during kickstage burn is large but remains nearly constant until the absolute value of the 
desired DLA exceeds the SS orbital inclination. A major portion of this velocity loss is 
caused by the trajectory constraints . The outbound Tv^ burn puts the kickstage and 
payload on an escape hyperbola. During this burn and the long coast phase that follows, 
the altitude increases very rapidly and the kickstage burn occurs at high altitudes, re- 
sulting in large losses. For DLA's whose absolute value is larger than the SS orbital 
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inclination, some inclination change is made during the kickstage burn and the losses 
increase rapidly. It should be pointed out that the velocity added to change orbital in- 
clination in order to meet mission DLA requirements is not a loss. However, because 
of the definition of ideal impulsive mission velocity (which is independent of the DLA) , 
the velocity required to change orbital inclination is included in the velocity losses along 
with the loss due to gravity and other trajectory -shaping losses. Outbound and inbound 
losses vary with DLA; however, the sum of these losses remains nearly constant. The 
trade-off between outbound and inbound velocity losses for the two solutions is obvious 
from this figure . 

2 2 

The same data discussed in the preceding paragraphs for an energy of 50 km /sec 

2 2 

are presented in figures 8 to 13 for an energy of 100 km /sec . Behavior of the data 
is very similar and will not be discussed in detail. Again the payload loss for a range 
of DLA's between -28.5° and +28.5° is small, less than 1 percent of payload. 

The data presented in figures 14 to 17 were derived for an outer planet mission. 
They show the effect on payload of variations in selected mission. Tug, and kickstage 
parameters. Since the payload curve is symmetric about 0° DLA, the data presented 
are restricted to the solution giving maximum payload and positive DLA's. Data are 
provided for DLA's of 0°, 10°, 20°, and 30°. 

Figure 14 shows the variation of payload with SS orbital inclination. All payload 
curves have a minimum at an initial orbital inclination of about 45°. Since at an in- 
clination of 90° the nodal shift is zero, the optimum flight profile is planar. Payload 
loss due to nodal shift then is the difference between the payload at a given inclination 
and the payload for the 90° inclination planar case. Maximum payload loss due to nodal 
shift for the sample mission (28.5° orbital inclination, 30° DLA, and 100-km^/sec^ 
mission energy) is 14 kilograms, or about 2. 2 percent of the payload. Also, since the 
nodal correction for the 90° initial SS orbital inclination is zero, the Tug payload capa- 
bility is maximum for this case. However, the SS is limited by range-safety launch 
azimuth constraints to a maximum orbital inclination of 56° for launch from the Eastern 
Test Range. The results presented in figure 14 were derived by assuming that the SS 
payload capability is constant and independent of orbital inclination. In the actual case, 
SS payload capability decreases as orbital inclination is increased, which would tend to 
reduce the indicated Tug/kickstage payload. Therefore, for initial SS orbital inclina- 
tions between 28. 5° and 56°, maximum performance (for DLA's between -28° and +28°) 
is obtained by departing from a 28. 5° inclined orbit. 

Data presented in figure 15 were obtained by varying the Tug/kickstage separation 
coast time. Payload capability is quite sensitive to this coast time. As coast time is 
increased from 2 minutes to 14 minutes, payload decreases from about 850 kilograms 
to 610 kilograms. This payload loss is caused by increased gravitational losses as 
coast time increases. Theoretically, optimal location for any burn is near perigee. At 
separation the Tug and kickstage are on an escape hyperbola. As coast time increases. 
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the two burns (retroburn and kickstage burn) will occur further and further from perigee. 
This, of course, increases gravitational losses and decreases payload. 

Another important mission parameter that affects payload is total trip time . Total 

trip time is defined as the time from Tug departure to return to the SS orbit. For short 

trip times the nodal correction is small, but the losses during Tug retroburn become 

large; for long trip times the nodal correction becomes large and the losses during 

retroburn become small. Therefore, there is an optimal trip time that gives maximum 

payload. Payload as a function of trip time is presented in figure 16. Optimum trip 

2 2 

time for this energy (100 km /sec ) is approximately 26 hours. However, the payload 
penalty for trip times between 20 and 36 hours is quite small (less than 1 percent of the 
payload) . 

Another Tug parameter investigated is the Tug thrust level, and the effect on pay- 
load is presented in figure 17. As initial thrust increases, the payload capability also 
increases because of decreased gravitational losses. These data are obtained by assum- 
ing a constant Tug engine weight. In the actual case, engine weight varies with thrust, 
which will affect the indicated payload capability . To select an optimum thrust level for 
the Tug, the variation of engine weight with thrust level would have to be included in 
determining payload capability . 

Kickstage burn (similar to the Tug retroburn) occurs on an escape hyperbola at high 
attitudes, and gravitational losses are large. These losses can be reduced by decreas - 
ing the kickstage burn time, which is equivalent to increasing the engine thrust level. 

The increase in payload capability that results from an increase in kickstage thrust is 
shown in figure 18. The results were obtained by assuming that the kickstage engine 
weight remains constant as its thrust varies. (The same assumption was made about 
Tug engine weight in determining the effect of Tug thrust on payload, as discussed pre- 
viously.) In this case, however, there is a one-to-one trade-off between engine weight 
and payload weight, and actual payload capability can be computed by subtracting the 
change in engine weight for a given thrust from the indicated payload capability. Opti- 
mum kickstage thrust level could be determined from this type of figure for a given mis- 
sion by including the effect of engine weight on payload. 


SUMMARY OF RESULTS 

In this report, mathematical equations were developed to obtain the maximum pay- 
load capability for a two-stage. Space Shuttle (SS) upper -stage configuration. This con- 
figuration consists of a reusable Tug and an expendable kickstage. The analysis 
assumes that correction for SS-orbit nodal precession is done by the Tug during the 
outbound burn and the retroburn. Planetary injection conditions are satisfied by a given 
energy and a given declination of the outgoing asymptote. The procedure presented may 
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be used to determine SS launch time to satisfy specific right ascension requirements for 
a given mission. The multipoint boundary value problem derived is solved by using a 
Newton-Raphson iterator. Partial derivatives required by this iteration technique are 
approximated by means of finite -difference techniques. 

To demonstrate the capability of the analytic technique, these equations were ap- 
plied to a vehicle configuration consisting of a cryogenic reusable Tug and a space- 
storable, pressure-fed, expendable kickstage. 

Payload as a function of the declination of the outgoing asymptote was computed for 

2 2 

mission energies of 50 and 100 km /sec and a set of nominal mission constraints. For 
these nominal mission constraints the maximum payload capabilities for the two mission 
energies were 2400 and 680 kilograms, respectively. This represents payload losses 
of 820 and 520 kilograms, respectively, when compared with the payload capability de- 
termined by using ideal impulsive computations without nodal correction. Therefore, 
it is obvious that ideal impulsive payload computations are inadequate and that the tra- 
jectories must be integrated. 

Perturbation analysis done for the outer -planet mission (vis-viva energy of 100 km^/ 
sec ) shows that payload is quite sensitive to Tug thrust as well as to kickstage thrust 
and Tug/kickstage separation coast time. Optimum total trip time is approximately 
26 hours; however, the payload loss for total trip times between 20 and 36 hours is 
quite small. Examination of departure -orbit inclination shows that maximum payload 
for an Eastern Test Range SS launch is obtained by departing from a 28. 5° inclined 
orbit. Of course, when high declinations are required, the departure -orbit inclination 
must be increased. 

The payload loss to correct for nodal precession is small (14 kg) for the optimal 
trajectories derived. The results further show that the energy range of a reusable Tug 
can be greatly extended by using this two-stage configuration and that significant payload 
capability is available for comet and asteroid types of missions and for some outer 
planet missions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 27, 1976, 

491-02. 
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APPENDIX A 


SYMBOLS 


a time at beginning of a phase, sec 

b time at end of a phase, sec 

C defined by eq. (18), kg-sec 

2 2 

E energy, m /sec 

e eccentricity 

f unit thrust direction 

3 2 

G gravitational constant of Earth, m /sec 

H Hamiltonian, kg 

H* augmented Hamiltonian, kg 

h angular momentum, m /sec 

h^ unit angular momentum vector of the return orbit 
1 inclination, rad 

-3 

J oblateness parameter (1.624x10 ) 

K kappa function defined by eq. (13), sec 

I right-handed Cartesian coordinate directions (Zj^ points to initial SS orbital line 
of nodes and I ^ to North Pole) 

m mass, kg-secV^i 

mpk kickstage propellant loading, kg 

p semilatus rectum, m 

q constraints 

Rg equatorial radius of Earth (6.37816x10® m) 
r radius, m 

T return time, sec 

Tp trip time, sec 

t time, sec 

At TugAickstage separation coast time, sec 

V velocity, m/sec 
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Vg exhaust velocity, m/sec 

V 1 Tug exhaust velocity, m/sec 

V „ kickstage exhaust velocity , m/sec 

e, z 

direction of outgoing asymptote 
AV„ apogee delta velocity, m/sec 

AV perigee delta velocity, m/sec 

AVj total delta velocity, AVrj. = AV^ + ^Vp, m/sec 

u argument of latitude, rad 

x,y dummy variables 

(X angle defined by eq. (B3.4), rad 

^ mass flow rate, kg/sec 

y adjoint variable associated with thrust direction, kg 

5 yaw attitude, rad 

e arbitrary constant 

^ dummy variable, sec 

7] true anomaly, rad 

6 defined by eq. (B3.3) 

X adjoint variable associated with velocity, kg -sec /m 

jU adjoint variable associated with position, kg-sec/m 

a adjoint multiplier associated with mass, sec 
T adjoint variable associated with state variable ^ , kg 

<p declination of outgoing asymptote, rad 

(// pitch attitude, rad 

longitude of ascending node, rad 
AS2 nodal precession 

w argument of pericenter , rad 

Ao; apsidal motion 

Subscripts: 

D desired 
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H hardware 

i index 

Superscripts: 

~ vector quantity 

unit vector quantity 
time derivative, d/dt 
Other notations: 

scalar product of two vector quantities 
X vector product of two vector quantities 

V- gradient of a scalar with respect to vector 

( I magnitude 

A difference 


APPENDIX B 


DERIVATION OF EQUATION (21.2) 

The gradients of declination with respect to radius and velocity have been derived 
in reference 1. They are given by 


Vjrto = { " [" a r + e] • i. 


+ 






a 


(^3 ■ r)r 


+ tan Of • 


e) 

h r - 


V 

h 




(Ig • ?)V + a ■ V)^3 + i (^3 • ^)VphJ|j (Bl) 


V^(p = I- [- tan O' r + &] - t. 
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Vye 


+ tan o -I — 
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?3 -(?3 • f)f 


VyhP “ 


h 


cos (p 


(B2) 


where 

sin (?? = i 3 ■ (B3 . 1) 

= cos Q? r + sin a & (B3 . 2) 

§ = h X r (B3 . 3) 

and 

a=cos"^-l-7? (B3.4) 

e 
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Adding the cross product of equation (Dl) and V to the cross product of equation (D2) 
and r gives 


X V + ^ ^ Of r + e] • i ||yjrn x r + VyT/ X vj 


+ [v^e X r + Vye x vj + x r - tan a — [(Zg • r)V x r + (r • V)f g x ^ 

e 


+ tan a - Tf , X V - (Z% • r)r x v1 - tan a — — - fv-h x r + V^h x vl|| .cos_a 

h L 3 ' 3 ^ J L r V cos cp 


(B4) 


It was shown in reference 6 , that VyX x V + V— X x r = 0 if X = X[r, V, r • In 
equation (B4) 77 , e, and h are of this form; therefore, the corresponding cross product 
sums are zero and equation (B4) reduces to 

X F + Vy<p X V = ^Zg X |tan a [V - (f • V)r] + rjj 

where the term rV/h [v - (r • V)r] = 0 and 

V-^xr + Vy<pxV = -J— ZgXV^ (B5) 

cos (p 
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TABLE I - NOMINAL MISSION CHARACTERISTICS 


2 2 

Range of vis-viva energy, km /sec 35 to 120 

Range of declinations, deg -34 to +34 

Total trip time, days ... . 1 

TugAickstage separation coast phase, min ... 10 

Circular departure -orbit altitude, km 185 

Departure -orbit inclination, deg 28.5 


TABLE n, - BASELINE TUG CHARACTERISTICS 


Total vehicle weight, kg 28 622 

Weight at Tug burnout, kg .' 2767 

Engine specific impulse, sec 460 

Engine thrust, N . . . 66 723 

Performance reserve, percent of specific impulse . . . . 2 

Kickstage adapter weight, kg ^ 45 


TABLE m. - KICKSTAGE CHARACTERISTICS 


Propellant weight, kg 2306 

Specific impulse, sec 378 

Stage burnout weight, kg 552 

Payload adapter weight, kg 9 

Engine thrust, N 13 345 






0° declination. 


Inclination at Tug kickstage Departure argument of latitude, deg 

separation, deg 



Declination. DLA. deg 


Figure 3. - Optimum Tug departure argument of latitude measured from Space Shuttle orbit 
ascending node for a mission energy of 50 km^/sec\ Initial SS orbital inclination, 28. 5°. 



Decimation. DLA. deg 


Figure 4 - i^perbolic orbit inclination at Tug/kickstage separation for a mission energy of 
50 km^/secs Initial Space Shuttle orbital inclination. 28.^. 



Pseudo right ascension, deg Velocity loss, m/sec 



Figure 7 - Velocity loss compared with ideal impulsive orbital 
an initial Space Shuttle orbital inclination of 28. ^ and a mii 
sec^ 



Figure 8. - Optimum pseudo right ascension of the outgoing 
rlnrkwkp frnm Soace Shuttle orbit ascend ina node for a m 


as 

is: 



nclination a 






Payload, kg Nodal change during outbound Tug burn, deg 




Decimation. DLA, deg 


Figure 13. - Velocity loss compared with ideal impulsive ortital transfer for departure 
from an initial Space Shuttle orbital inclination of 28 5° and a mission energy of 
lOOkm^/sec^. 


Declination, 






Figure 13. - Effect of Tug/kickstage separation coast time on pay- 
load. Initial Space Shuttle orbital Inclination. 28.^; mission 
energy, lOOkm^/sec^ 


Declination. 



Figure 16. - Effect of Tug total trq) time on payload. 
Initial Space Shuttle orbital inclination, 28.^; 
mission energy, 100km^isec‘. 
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Figure 17. - Effect of Tug initial thrust on payload - fixed Tug engine weight Initial 
Space Shuttle orbital inclination, 28.5°; mission energy, 100 km^/ sec’ 



J ^ I I ^ ^ I ^ I ^ I , 

4 8 12 16 20 24 28 32 36 40 44x10^ 

KIckstage thrust. IM 


Figure 18. - Effect of kickstage thrust oi^ayload - fixed kickstage engine weight. Initial 
Space Shuttle orbital Inclination, 28. mission energy. lOOkm^fsec^. 
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